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Introduction 


The  Eastern  Pacific  Hurricane  Strike  Probability 
Program  (EPSTRKP)  is  nearly  identical  in  concept  to  the  com¬ 
parable  Western  Pacific  Strike  Probability  Program  (STRIKP) 
(Jarrell,  1978)1.  There  are,  however,  some  differences  in 

detail.  An  earlier  study  on  eastern  Pacific  tropical  cyclones 

2 

was  conducted  by  W.J.  Thompson  and  R.L.  Elsberry  (1979)  . 
Thompson  and  Elsberry  subsequently  expanded  the  study  to  cover 
the  central  Pacific  (140°W  to  180°)  and  included  augmented 
tracks  of  those  cyclones  which  recurved  over  North  America. 
This  latter  step  was  considered  to  be  necessary  because  the 
earlier  geographical  distribution  of  errors  appeared  to  be 
unrealistically  low  along  the  coast.  This  was  attributed 
to  the  absence  of  both  forecast  positions  (since  overland 
dissipation  was  forecast)  and  verification  positions  (actual 
dissipation)  overland,  thus  only  correctly  forecast  overwater 
tracks  were  verifiable.  The  augmentation  of  recurvature 
tracks  provided  an  increase  in  verifying  positions,  but  made 
little  or  no  change  in  the  geographical  error  distribution. 

Original  Study 

Thompson  and  Elsberry* s  original  study  was  based 
on  the  official  forecasts  issued  by  the  Weather  Service  Fore¬ 
cast  Office,  San  Francisco  for  the  years  1971  to  1977.  It 
included  only  those  forecasts  in  the  San  Francisco  area  of 
responsibility,  which  represents  over  80%  of  the  eastern 


Darrell,  Jerry  D. ,  Tropical  Cyclone  Strike  Probability  Fore¬ 
casting.  NAVENVPREDRSCHFAC  Contractor  Report  CR78-01. 
December  1978. 

2 

Thompson,  W.J.,  and  R.L.  Elsberry,  A  statistical  analysis  of 
eastern  Pacific  tropical  cyclone  forecast  errors.  Twelfth 
Technical  Con.  on  Hurricanes  and  Tropical  Meteorology,  New 
Orleans,  April  1979. 


Pacific  tropical  cyclones.  After  the  forecasts  for  the  year 
1978  were  added  there  were  2036  verifiable  24-hour  forecasts. 

Thompson  and  Elsberry  also  performed  discriminant 
analyses  on  forecast  error.  Some  of  the  results  of  those 
analyses  will  be  presented  in  Appendix  A. 

Augmented  Study 

Thompson  and  Elsberry' s  1971-78  data  base  was  aug¬ 
mented  by  extending  the  postanalysis  tracks  of  hurricanes  over 
North  America  as  far  as  possible.  This  was  done  with  the  aid 
of  twice  per  day  surface  weather  charts  and  satellite  mosaics. 
The  former  were  of  limited  use  because  the  rough  terrain  of 
western  Mexico  makes  interpretation  of  surface  reports  in 
the  vicinity  of  a  hurricane  exceedingly  difficult.  The 
satellite  mosaics  were  the  roost  usable  information  source 
except,  of  course,  there  is  little  precision  in  locating  a 
cyclone  center  when  it  has  degenerated  into  a  cloud  mass  over 
rough  terrain.  Nevertheless  the  tracks  of  those  that  were 
extended  were  reasonably  certain. 

The  original  data  base  was  also  augmented  by  adding 
those  forecasts  for  the  central  Pacific  issued  by  the  Central 
Pacific  Hurricane  Center  located  at  the  Weather  Service  Fore¬ 
cast  Office,  Honolulu. 

Statistics  for  the  study  before  and  after  track 
extension  and  central  Pacific  augmentation  are  shown  in 
Table  1. 


The  extension  of  overland  tracks  increased  the 
number  of  verifiable  forecasts  from  1%  at  24  hours  to  6%  by 
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72  hours.  This  had  the  net  effect  of  increasing  average 
forecast  errors  as  was  expected. 

Figure  1  shows  average  official  24-hour  movement 
forecast  errors  on  a  geographical  plot.  Figure  lb  shows  the 
errors  for  the  earlier  study  forecasts  while  figure  la  shows 
the  same  information  after  augmentation  by  extension  of  over¬ 
land  verification  tracks  and  the  inclusion  of  central  Pacific 
forecasts.  Except  for  the  obvious  westward  continuation  of 
contours,  the  two  figures  are  virtually  identical. 


24  Hour 

48  Hour 

72  Hour 

A 

B 

c 

A 

B 

C 

A 

B 

C 

S-N  Error 

Mean 

-7 

-8 

-5 

-12 

-16 

-10 

-15 

-25 

-16 

Std.  Dev. 

87 

87 

86 

154 

155 

154 

213 

219 

217 

W-E  Error 

Mean 

-6 

-7 

-8 

-  2 

-  9 

-15 

-  6 

-30 

-33 

Std. Dev. 

104 

104 

102 

197 

200 

194 

295 

310 

303 

Vector  Error 

Mean 

115 

115 

112 

215 

218 

212 

312 

326 

318 

Std.  Dev. 

74 

74 

73 

128 

129 

130 

188 

198 

198 

Corr  Coef 


S-N  vs  W-E 

.05 

.05 

.09 

.08  .08  .16 

.18  .18 

.24 

Cases 

2036 

2058 

2327 

1276  1309  1541 

924  976 

1163 

Table  1. 

Evolution 

of  error 

statistics  from  1971-78. 

Means  and 

standard  deviations  are  in  n  mi. 

A  Eastern  Pacific 

B  Eastern  Pacific  with  overland  track  extensions 
C  Eastern  and  central  Pacific  with  overland  track  extensions. 
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The  addition  of  central  Pacific  forecasts  increased 
the  number  of  verifiable  forecasts  by  13,  18  and  19%  at  24, 

48  and  72  hours  respectively.  The  average  vector  error  de¬ 
creased  with  the  addition  of  these  forecasts,  although  the 
bias  (component  mean  errors)  remained  about  the  same  in 
magnitude  but  its  orientation  became  more  westerly. 

Operational  Products 

Two  distinct  strike  probability  products  will  be 
available  under  operational  evaluation  during  the  1980 
hurricane  season  for  the  eastern  North  Pacific. 

Product  1 .  Tropical  cyclone  strike  probabilities 
for  preselected  points.  This  can  be 
disseminated  automatically  to  a  dis¬ 
tribution  list  by  Fleet  Numerical 
Oceanographic  Center  (FNOC)  via  AUTODIN 
initially  and  possibly  later  via  the 
Automated  Weather  Network  (AWN).  In¬ 
cluded  in  this  product  will  be  a  fore¬ 
cast  class  specification  for  confidence 
estimates  for  the  Naval  Western  Oceano¬ 
graphy  Center  (NWOC)  Pearl  Harbor 
(refer  to  figure  3  of  Appendix  A). 

Product  2.  Individual  user  requests  for  tropical 
cyclone  strike  probabilities  via  the 
Automated  Product  Request  (APR)  system 
(AUTODIN  only). 


Product  1  could  be  generated  routinely  by  FNOC 
upon  receipt  of  the  NWOC  Pearl  Harbor  tropical  cyclone  warning 
every  six  hours.  Product  1  would  give  the  probabilities  of  a 
particular  tropical  cyclone  being  within  75  n  mi  (left)  or 


50  n  mi  (right)  relative  to  forecast  track  of  nine  preselected 
points  of  interest.  Although  subject  to  change  the  points 
currently  listed  within  the  program  are:  Acapulco,  Mazatlan, 
Puerto  Vallarta,  La  Paz,  San  Diego,  Hilo,  Honolulu,  Johnston 
Island  and  Midway  Island.  The  strike  probabilities,  computed 
upon  receipt  of  each  6-hourly  warning  and  given  at  12-hour 
intervals  after  warning  time,  are  presented  in  two  forms.  The 
first  is  the  instantaneous  probability,  valid  at  a  single 
instant  of  time  only.  The  second  is  a  time  integrated 
probability  —  the  probability  that  a  strike  will  occur  at 
some  time  between  the  effective  time  of  the  warning  and  mul¬ 
tiples  of  12  hours  thereafter.  Similarly  probabilities  of 
30  and  50  kt  winds  are  expected  to  be  added  to  this  message 
at  a  later  date. 

Product  2  would  be  run  only  upon  request.  The  user 
would  make  his  request  to  FNOC  via  AUTODIN.  He  would  include 
information  sufficient  to  identify  the  tropical  cyclone,  the 
point  of  concern  ( 1 at it ude/ longitude ) ,  and  the  radii  about 
that  point  describing  the  area  considered  to  constitute  a 
strike.  The  output  would  be  in  the  same  form  as  in  product 
1  (i.e.,  instantaneous  and  time  integrated  strike  probabilities 
at  12-hour  intervals  after  warning’  time) .  *•  ....... 

An  example  follows  to  show  the  user  how  the  output 
will  appear.  The  example  is  Tropical  Storm  ANDRES  at  1800 
GMT  2  June  1979.  At  this  time  ANDRES  was  175  n  mi  south  of 
Acapulco,  Mexico  with  35  kt  winds.  It  was  expected  to  move 
northwest  at  8  kts  for  the  first  12  hours  becoming  westerly 
and  finally  westsouthwest  during  June  5th  (GMT).  Its  intensity 
was  expected  to  increase  over  the  3-day  period  to  55  kts. 


-7- 


Two  Eastern  Pacific  Strike  Probability  Program 
(EPSTRKP)  runs  for  ANDRES  are  discussed  below. 

Run  1  is  a  PNOC  originated  run  (Product  1 )  at 
02 /1800  GMT. 

Run  2  is  in  response  to  a  hypothetical  user  also 
at  02/1800  GMT  specifying  an  area  within  100  n  mi  of  a  point 
(15°N,  102°W).  His  request  would  have  gone  to  FNOC  via 
AUTODIN  message  as  an  APR  formatted  message  (see  Table  2). 
Required  input  is  at  least  one  Area  of  Concern  (lat/long) 
and  radii  to  the  left  and  right  of  that  point  (relative  to 
forecast  motion). 

Tables  3  and  4  illustrate  the  output  from  Runs  1 
and  2,  respectively.  These  tables  also  contain  some  descrip¬ 
tive  information. 


BT 

UNCLAS / / N03 1 60 / / 

TROPICAL  CYCLONE  STRIKE  PROBABILITY  REQUEST,  EASTERN  PACIFIC 
Q92X0001 

/APR, AP(EPSTRKP) , (other  entries  on  this  line  as  required) 
/STM , NM( ANDRES ) , NR( EP02 ) ,  DH( 7906021800 ) / 

/ AOC , LA( 1 50N ) , LO( 1020W ) , RL( 100 ) , RR( 100 ) / 


(as  many  AOC  lines  as  needed) 


/AAD, 

etc.  (as  needed) 

/PARA , 

/ERK/  (required  end) 

BT 

/STM:  Storm  line 

NM:  Name  of  cyclone 

NR:  Cyclone  number,  Ocean  Basin  EP 

DH:  Effective  Date/time  of  warning. 

(Day  02  hour  1800  GMT) 

,/AOC:  Area  of  concern  line 

LA:  Latitude  of  point  of  concern.  LA(150N)=15.0°  north. 

LO:  Longitude  of  point  of  concern.  LO(1020W)=102.0°  west. 

RL:  Radius  of  area  of  concern  to  left  of  storms  track. 

RR:  Radius  of  area  of  concern  to  right  of  storms  track. 

Usually  RL  is  greater  than  RR.  Default  values  of  75/50  nm 
will  be  used  if  both  RL  and  RR  are  zero  or  blank. 


Note:  One  input  record  will  be  written  for  each  /AOC  (including 
storm  information).  Request  message  in  accordance  with 
FLENUMOCNCEN ,  1977:  ASWOCAS  Request  Procedures  Manual, 

Vol .  2. 

Table  2.  Sample  Automated  Product  Request  (APR)  System  Message. 


=  Eastern  Pacific 

DH( 7906021800)  = 
021800Z  June  1979 
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Run  1  Output  (Product  1) 


STRIKE  PROBABILITY  FORECASTS 
ANDRES  021800Z 

ACAPULCO  00ININ*1 20815  240418  360218  480119  600119  720119 

MAZATLAN  00ININ  12ININ  24IN1N  36ININ  48ININ  60IN01  72IN02 

P  VALLRTA  00ININ  12ININ  24ININ  360101  480104  600106  720107 

LA  PAZ  00ININ  12ININ  24ININ  36ININ  48ININ  60ININ  72IN01 

SAN  DIEGO  THREAT  NIL* 

HILO  THREAT  NIL 

HONOLULU  THREAT  NIL 

JOHNSTN  I  THREAT  NIL 

MIDWAY  I  THREAT  NIL 

FOR  NWOC  PEARL . . CLASS  =  THREE 

PROBABILITIES  BASED  ON  FOLLOWING  FORECAST 


FORECASTS:  Time  12  hr  Latitude  14. 8N  Longitude  100. 4W  Max  Wind  40  kt 


LAT/LONG  of  preselected  points  are  stored  within  program. 
Strike  is  predefined  to  occur  if  tropical  cyclone  passes 
within  75  n  mi  radius  (left)  or  50  n  mi  radius  (right)  of 
track  of  tropical  cyclone. 


♦THREAT  NIL  means  all  probabilities  for  this  station  were 
<0.5%.  IN  means  insignificant  (<0.5%). 


Table  3.  Output  from  Run  (1). 
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Run  2  Output  (Product  2) 

STRIKE  PROBABILITIES  FOR  TROPICAL  CYCLONE  ANDRES 
FROM  021 800 Z  BASED  ON  FOLLOWING  FORECAST 

001400989035  121481004040  241501020045  481421049050  721321070055 
STRIKE  IS  BEING  WITHIN  100NM  RIGHT  AND  100NM  LEFT  OF  15. ON  102. OW 
STRIKE  PROBS  000202  123840  243143  361443  480743  600443  720343 

If  TIME 

/  PR0B(%)  that  ANDRES  will  be  in 

/  area  at  051800  (Warning  time  + 

/  72  hours)  was  3% 

/  PROB  that  ANDRES  will  be  in  - 

/  area  some  time  between 

/  021800Z  and  051800Z  (72  hour 

/  period)  was  43%. 

ABBREVIATIONS : 

Number  01-99;  strike  probability  in  % 

IN  =  insignificant;  p<0.5%  Prevents  representation  of  0% 

and  100%  which  occur  only  as  an 
approximation . 

The  input  forecast  data  is  error  checked  only  in  that  the  tropical 
cyclone  forecast  motion  is  computed  between  forecast  points.  If  vector 
motion  deviates  substantially  from  the  climatological  mean,  the  following 
warning  message  will  appear  in  all  products: 

***  UNUSUAL  MOTION  —  PLEASE  RECHECK  WARNING  DATA  *** 


TABLE  4.  OUTPUT  FROM  RUN  (2). 
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APPENDIX  A 


Discriminant  Analysis 

The  discriminant  analyses  routine  of  the  UCLA  BIOMED 
(Dixon,  1975)  series  was  run  on  the  data  to  develop  functions  to 
discriminate  on  forecast  error.  Predictands  were  forecast  error 
group  numbers  1,  2  and  3.  The  groups  were  determined  by  using 
in  turn  the  24-,  48-  and  72-hour  forecast  errors  to  split  the 
forecasts  into  three  equal  groups  according  to  error  magnitude; 
group  1  had  the  smallest  errors  and  group  3  the  largest.  Table  A- 
1  shows  average  forecast  errors  for  the  three  error  classes, 
where  the  classes  were  imperfectly  discriminated  by  applying 
functions  developed  on  each  of  the  predictands.  If  discrimina¬ 
tion  were  perfect,  the  groups  and  classes  would  be  identical. 


Applied  to:  24  hr  errors  48  hr  errors  72  hr  errors 


Class: 

Used  for  Split: 

i 

2 

3 

1 

2 

3 

1 

2 

3 

24-hour  errors 

88 

107 

138 

185 

214 

252 

280 

335 

371 

48-hour  errors 

95 

98 

132 

174 

194 

265 

266 

301 

402 

72-hour  errors 

95 

94 

128 

176 

187 

255 

265 

288 

393 

Table  A-l.  Average  movement  forecast  errors  for  24,  48  and  72  hours  classed 
by  discriminant  functions.  The  functions  were  derived  on  fore¬ 
casts  of  one  time  length,  but  applied  to  the  other  lengths. 


Based  on  the  difference  in  group  means,  the  discrimination 
provided  by  the  functions  developed  on  48  hour  errors  (under¬ 
lined)  is  superior  to  the  other  two.  There  is  generally  poor 
discrimination  between  classes  1  and  2  as  evidenced  by  the  close¬ 
ness  of  their  means.  Class  three,  the  difficult  forecast, 
appears  to  be  well  separated  from  the  other  two  classes. 

Eight  predictors  were  selected  from  23  candidates.  Table  A- 
2  defines  the  eight  predictors,  gives  their  means,  standard 

^Dixon,  W.  J.,  BMDP  biomedical  computer  programs.  University  of  California 
Press,  Berkeley,  1975. 
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deviations  and  the  coefficients  that  define  the  discriminating 
functions.  A  close  examination  of  Table  A-2  suggests  less 
than  optimum  conditions  for  discriminant  analysis.  Some  of 
the  predictors  are  clearly  not  normally  distributed,  (i.e., 
direction  is  strongly  bimodal),  some  are  closely  related 
(direction  and  adjusted  direction).  These  are  not  serious 
since  the  worst  offender,  direction,  is  one  of  the  least 
significant  predictors  (see  diagram  at  the  origin  of  figure  A- 
1  where  length  of  the  vector  relates  to  predictor  contribution 
to  error  discrimination).  Other  aspects  of  the  discriminant 
analysis  that  might  appear  questionable  are  the  small  dif¬ 
ferences  in  predictor  means.  Sample  size  was  2327  cases 

so  each  class  has  on  the  order  of  seven  hundred  cases. 

Differences  in  means  of  1 . 96St . Dev . / /  to  are  significant  at  the  5% 
level,  which  means  differences  greater  than  0.07S  are  sig¬ 
nificant.  For  example,  the  standard  deviation  for  variable 
#8  is  about  1.5  kts ,  therefore  differences  in  the  mean  of  .11 
kts  would  be  significant.  The  actual  differences  are  .12  kts 
and  .51  kts  between  the  class  1  mean  and  the  other  two  and 
.39  kts  between  classes  2  and  3,  all  are  significant.  Of 
course  not  all  the  differences  in  the  means  are  significant 
at  the  5%  level,  but  most  are. 

One  comparable  way  to  look  at  the  differences  in 
standard  deviations  is  to  perform  the  same  test  on  the  means 
as  above,  but  use  each  of  the  class  standard  deviations. 
Intuitively,  if  there  is  no  difference  in  the  outcome,  the 
difference  in  the  standard  deviations  is  not  important.  There 
are  very  few  comparisons  where  the  outcome  depends  on  which 
of  the  three  standard  deviations  is  used. 

Figure  A-l  shows  an  x-y  space  with  function  1  on  the 
x-axis  and  function  2  the  y-axis.  The  values  of  these  functions 
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Figure  A-l.  Function  1  vs  Function  2  in  x-y  space.  +'s  show  location 
o£  intersection  of  functions  evaluated  at  the  three  class  means.  The 
space  is  partitioned  according  to  class.  The  diagram  at  the  origin 
depicts  the  effect  of  an  increase  of  one  standard  deviation  in  each 
of  the  eight  predictors.  Star  refers  to  Tropical  Storm  ANDRES  on  2 
June  1979  at  1800  GMT.  (See  sample  output,  Table  3  of  text.) 


have  been  evaluated  at  the  group  mean  for  each  predictor  (see 
Table  A-2).  These  mean  points  are  plotted,  and  the  space  is 
divided  into  three  areas  whose  member  points  would  be  nearest 
to  the  included  mean  point.  Since  there  is  little  difference 
in  class  1  and  2,  the  far  left  of  the  diagram  is  the  difficult 
area,  with  the  right  and  center  the  average  forecasts.  The 
origin  represents  the  two  functions  evaluated  at  the  population 
mean,  versus  class  mean,  of  each  predictor. 

The  diagram  at  the  origin  represents  the  effect  of 
one  predictor  alone  being  increased  by  one  standard  deviation. 
This  provides  some  insight  into  the  causes  of  forecast  dif¬ 
ficulty  in  EASTPAC.  See  for  example  the  effect  of  a  deviation 
of  predictor  5,  Fix  Accuracy.  This  is  actually  the  mean  warning 
position  error  for  all  warnings  based  on  combinations  of 
methods  used  to  locate  the  center  (aircraft,  satellite,  radar, 
etc.).  While  the  mean  error  is  around  25  n  mi,  the  worst 
method,  extrapolation  alone,  has  a  50  n  mi  average  error, 
which  is  about  3  standard  deviations  away  from  the  mean. 

Notice  in  figure  A-l,  three  lengths  of  the  small  vector  labeled 
5  moves  a  forecast  from  the  average  into  the  difficult  (class 
three  area)  forecast  region  by  a  large  margin. 

The  star  in  figure  A-l  represents  an  actual  computation 
of  Function  1  (x=-2.47)  and  Function  2  (y=0.60)  for  tropical 
storm  ANDRES  on  2  June  1979.  (ANDRES  was  discussed  as  the 
output  example  in  the  operational  products  section  of  this 
report).  We  expect  this  is  a  difficult  forecast  and  will 
likely  have  a  large  error.  The  discriminant  analysis  predic¬ 
tors  and  basic  variable  are  given  below  as  are  the  terms 
of  Function  1  (x)  and  Function  2  (y). 


Basic  Variable 


Predictor 


Function  1 


Function  2 


1 

Direction 

NW=315° 

315° 

-1.58 

2.21 

2 

Latitude 

14°N 

14  # 

-1.46 

1.47 

3 

Max  wind 

35  kt 

35  kt 

0.53 

-0.74 

4 

Adj  dir 

315-110 

205° 

-3.08 

2.46 

5 

Fix  acy 

Sat /Ext rap 

43.5  n  mi 

-2.52 

-0.35 

6 

Jul  date 

153 

54.76 

-1.10 

-1.15 

7 

Adj  Long 

98. 9°W 

20.10° 

0.86 

1.41 

8 

North  mvmt 

5.66  kts 

2.66  kt 

-0.84 

0.24 

Constant 

6.72 

-4.95 

-2.47 

0.60 

Generally  the  terms  which  are  contributing  to  the 
forecast  difficulty  are  those  where  Function  1  is  negative 
(direction,  both  terms,  and  fix  accuracy)  and  to  a  lesser 
extent  those  negative  terms  in  Function  2  which  are  not  of 
consequence  here. 

Figure  A-2  shows  plots  of  24-,  48-  and  72-hour  unit 
probability  ellipses.  A  point  on  a  unit  ellipse  is  one  stan¬ 
dard  deviation  away  (+)  from  the  mean  in  one  coordinate  when 
the  other  is  at  the  mean.  A  unit  probability  ellipse  is 
equivalent  to  a  39%  probability  ellipse. 
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Figure  A-2.  Unit  probability  ellipses  of  forecasting  errors  from  three  forecast 
difficulty  classes  at  24,  48  and  72  hours.  The  class  parameters 
are  given  in  the  inset  table.  The  origin  is  the  forecast  position. 
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The  three  nested  sets  represent  the  three  classes 
of  errors.  The  nested  ellipses  represent  in  increasing  size 
the  24-,  48-  and  72-hour  forecast.  The  class  2  ellipses  are 
larger  than  those  of  class  1  but  the  difference  is  small.  In 
general  the  class  3  ellipses  are  similar  in  shape  and  orien¬ 
tation  to  those  of  classes  1  and  2  but  corresponding  ellipses 
are  roughly  twice  as  large  (in  area).  The  inset  table 
summarizes  the  ellipse  parameters. 

The  discrimination  can  be  contrasted  to  that  reported 
2 

by  Nicklin  (1977)  for  the  western  Pacific.  In  that  case 
the  24-hour  errors  were  far  better  discriminated,  but  the 
72-hour  errors  were  not  discriminated  as  well.  It  can  be  seen 
from  Table  A-l  that  if  the  24-hour  error  had  been  used  as  the 
basis  for  discrimination  (as  it  was  in  the  western  Pacific) 
the  effect  would  have  been  to  make  the  results  more  similar. 

The  48-hour  discrimination  is  about  the  same  in  both  ocean 
basins.  One  difference  is  that  in  the  western  Pacific  the 
less  difficult  third  of  the  forecasts  were  well  separated 
from  the  remainder  whereas  in  the  eastern  Pacific  the  most 
difficult  third  of  the  forecasts  are  more  easily  isolated. 
Overall,  the  discrimination  appears  slightly  better  in  the 
eastern  Pacific,  perhaps  because  of  the  inclusion  of  infor¬ 
mation  relative  to  initial  position  accuracy  (see  variable 
5,  Table  A-2). 

Forecast  confidence  estimates  can  be  inferred  for 
each  forecast  class.  These  can  either  be  expressed  as  a 
percentage  of  occasions  when  the  actual  forecast  error  will 

Nicklin,  Donald  S. ,  A  Statistical  Analysis  of  Western  Pacific 
Tropical  Cyclone  Forecast  Errors,  M.S.  Thesis,  U.S.  Naval 
Postgraduate  School,  Monterey,  CA,  June  1977. 
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lie  between  zero  and  some  set  distance  (i.e.,  the  probability 
the  error  is  less  than  100  n  mi)  or  as  a  distance  which  will 
exceed  the  actual  error  with  a  set  probability  (the  radius 
of  a  75%  probability  circle).  Figure  A-3  presents  this  prob¬ 
ability  information  for  24-,  48-  and  72-hour  forecasts  (curves 
left  to  right)  and  for  classes  1,  2  and  3  (top  to  bottom). 

For  example  200  n  mi  represents  the  radius  of  95%,  57%  and  30% 
probability  circles  for  class  1  24-,  48-  and  72-hour  forecasts 
respectively.  Similarly  the  80%  confidence  limit  on  class  2 
forecast  errors  is  153,  296  and  465  n  mi  at  24,  48  and  72 
hours  respectively. 
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